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ABSTRACT 

We present results from our 830 ksec observation of the globular cluster 47 Tu- 
canae with the Chandra X-ray Observatory's High Resolution Camera-S. We 
limit our analysis here to the 19 previously known, localized millisecond pulsars 
(MSPs) in the cluster. This work more than doubles the sample of X-ray-detected 
MSPs observed with sensitivity to rotational variability; it is also the first survey 
of a large group of radio-discovered MSPs for which no previous X-ray pulsations 
have been detected and is therefore an unbiased survey of the X-ray properties 
of radio-discovered MSPs. We find that only 47 Tuc D, O and R show signifi- 
cant pulsations at the ^ 4-cr level, but there is statistical evidence for rotational 
variability in five additional MSPs. Furthermore, we constrain the pulsed mag- 
netospheric emission of 7 more MSPs using Monte Carlo simulations. The result 
is that the majority of the 47 Tuc MSPs are characterized by low pulsed frac- 
tions, ^ 50%. In cases where larger pulsed fractions are measured the folded 
pulse profiles show relatively large duty cycles. When considered with previous 
spectroscopic studies, this suggests that the X-ray emission arises from the neu- 
tron star's heated polar caps, and in some cases, from intra-binary shocks, but 
generally not directly from the star's magnetosphere. We discuss the impact of 
these results on our understanding of high energy emission from MSPs. 
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Subject headings: globular clusters: individual (47 Tucanae) — pulsars: general 
— stars: neutron — X-rays: stars 



1. Introduction 

Millisecond pulsars (MSPs) are old neutron stars spun-up by accretion of m ass and 



angular momentum from the matter of a donor binary companion (lAlpar et al.lll982l ). When 
compared to the canonical radio pulsar population they are distinguished by short spin 
periods, P ^ 25 ms, small spin-down rates, P ^ lCT 20 s/s, and thus low inferred dipole 
magnetic field strengths, Pdipoie oc (PP) 1 / 2 ~ 10 8 ~ 10 G, with large characteristic ages, r = 
P/2P ^ 1 Gyr. Studies of the ~ 150 known MSPs are difficult at wavebands outside of 
the radio due to their intrinsic faintness. The vast majority (~ 80%) of MSPs have binary 
companions that dominate at optical wavelengths, thus X-rays are an important avenue for 
studying MSPs. 

Currently, only 16 MSPs outside of 47 Tucanae (NGC 104, hereafter 47 Tuc) have been 
detected in X-rays, and only 12 of these have been observed with sufficient time resolution to 
explore variability on rotational timescales (see Tabled]). There are several proposed physical 
mechanisms capable of generating X-rays from these MSPs. Non-thermal emission processes 
in the neutron star magnetosphere generate power-law components in their X-ray spectra 
with characteristic photon indices T » 1.5-2. Pulsars in this class (e.g. PSRs B1937+21 and 
B1821— 24A; those above the first horizontal line in Table [1]) have large spin-down energies 
(E ^ 10 35 erg s -1 ), bright X-ray emission (L x ^ 10 32 erg s _1 ), low duty cycles and pulse 
profiles that closely resemble the radio emission in morphology and phase with large pulsed 
fractions, f p ^ 50% (see §3. lj) . Power-law spectral components can also be produced when 
the wind from the MSP interacts with material from the binary companion causing an intra- 
binary shock. These pulsars have similar properties to those above, with the exception 
that they lack strong rotational modulation (e.g. the 'black widow' pulsar, PSR 1957+20). 
Finally, heating of the neutron star polar cap by the bombardment of relativistic particles 
provides a mechanism for producing thermal X-ray emission. MSPs dominated by thermal 
spectra (e.g. PSRs J0437— 4715, J2124— 3358; most of those below the first horizontal line 
in Table [1]) are characterized by lower spin down energies (E <^ 10 34 erg s -1 ), lower X- 
ray luminosities (L x <^ 10 32 erg s -1 ) and pulse profiles with large duty cycles. As seen in 
Table [IJ the pulsed fractions of these MSPs are usually poorly constrained, but generally 
show f p « 50%. The emission of thermal cooling X-rays from the neutron star surface and 
those from pulsar wind nebulae are not thought to be important for these old objects, so we 
will not consider them here. 
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The unprecedented spatial resolution of Chandra has enabled detailed studies of MSPs 



vided the first census of low-luminosity X-ray sources in these clusters ( 


Rutledee et al. 


2004 


Becker et al. 2003: 


Bassa et al. 2004; Grindlav et al. 2002: Ransom et al. 


2004h. However, 


the largest endeavor has been Chandra 's observing campaign of 47 Tuc ( 


Grindlav et al. 


2001 



20021 ; iHeinke et al.ll2005t iBogdanov et al.ll2006l ). This work has shown that the spectral char- 
acteristics of the 47 Tuc MSPs are relatively homogeneous. Their luminosities fall in a narrow 
range, Lx ~ io 30-31 erg s _1 , and are well described by thermal spectral models with small 
emission radii, R e s ~ 0.1-3 km and temperatures of T e g ?5il-3xl0 6 K. The only exceptions 
are the radio-eclipsing binaries 47 Tuc J, O and W, which require additional power-law 
components above 2keV with T « 1.0-1.5. These results are reinforced by the findings of 
detailed spectroscopic studies by Chandra and XMM-Newton that have emphasized the dom- 
i nant thermal components in ne arby MSP X-ray spectra over the fainter power-law features 
(jZavlin et al.l |2002| ; IZavlinl 120061 ) . Consequently, we expect the predominant X-ray emission 
from the MSPs in 47 Tuc to arise from heated polar caps, and to be mo dulated just by 



rotat i ng a small area relative to the observer (i.e . sinusoidal pulse profiles; iGrindlay et al. 
20021 : Idieng & Taan]l2003l : IBogdanov et~aH l2006h . 



In this paper, we present high time resolution data capable of exploring the fast time 
variability of the X-ray counterparts to the 47 Tuc MSPs. In § [2] we present the details 
of the observations and data reduction, followed by the variability analysis (§ [3]) and an 
examination of the accuracy of the HRC-S time tags (§HJ). We find that the HRC-S remains 
capable of detecting fast modulation, but the MSPs in 47 Tuc lack strong variability on all 
timescales probed. We discuss the impact of this result on our understanding of the X-ray 
emission from MSPs in § [51 



Observations and Data Reduction 



Observations of 47 T uc were performed with Chandra's HRC-S detector (jZombeck et al. 



19951 ; iMurray et al.lll998l ). They began 2005 Dec 19 7:20 UT with 14 subsequent visits over 
the next 20 days for a total of 833.9 ks of exposure time (see Table [2] for a summary). The 
observing plan of dividing the 833.9 ksec into 50-100 ksec visits spread out over 20 days was 
adopted - instead of the optimal choice for pulsar detection of an uninterrupted observation - 
to mitigate a thermal limitation in space-craft operations. The HRC-S has a timing resolution 
of 15.625 ^s in the nominal energy range of 0.1-2keV, although it has essentially no energy 
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resolution. The data were analyzed using Chandra Interactive Analysis of Observations^ 
software (CIAO) version 3.3 and CALDB version 3.2.1. 

We began the data analysis of each observation by registering it to the first pointing 
(ObsID 5542) using the relative astrometry of the four brightest sources in the field. These 
corrections typically resulted in <^ Of! 5 corrections to the native astrometry. Data were 
filtered on pulse invariant (PI) channel to maximize the signal-to-noise ratio (SNR) for the 
known MSPs using the following approach. We compared the PI channel distribution of all 
counts extracted from regions corresponding to the known MSPs with that extracted from 
a background region. The background and source PI distributions are identical below PI = 
25 (within 1% in counts/area). Above this value the source region counts are significantly in 
excess, so we adopt PI = 25 as the lower PI limit for background filtering. We determined 
the upper-limit in PI channel by maximizing the total source SNR, finding the maximum 
when data at PI > 120 are excluded from analysis. Thus, we adopt a PI range of 25-120 
for all analyses; this decreases the MSP source count rates by 6%, while excluding 51% of 
the total background counts. The resulting background contribution is 17.9 counts arcsec -2 
over the entire observing span. 

Currently, positions are known for 19 of t he 22 MSPs in 47 Tuc. P ublished pulsar timing 



solutions can account for 16 of these positions (IFreire et al.ll200ll . l2003l ). while 47 Tuc R and Y 



have unpublished solutions (Freire et al. in preparation). 47 Tuc W has only a preliminary 
timing soluti on, but was localized by the eclipses of its optical counterpart (Freire et al. 



unpublished; Edmonds et al.ll2002l ). The close proximity of pulsars 47 Tuc G and I and 
47 Tuc F and S (separations 0'/12 and 0'.'7, respectively) do not allow them to be resolved 
by Chandra, so their counts must be considered jointly. For analysis purposes we attribute 
50% of detected photons to each pulsar (see §3.11) . 

We extracted photons from within circular regions surrounding each of the known MSP 
positions. The size of each extraction region can be found in Table [31 The radius was chosen 
adaptively in order to maximize SNR, but constrained to mitigate contamination by nearby 
objects. However, some contamination due to source crowding is unavoidable. We estimate 
this contribution by modeling the PSF as a Gaussian with o = 0.29", calculating the number 
of photons that fall in the extraction region of a given source from each neighboring source, 
and using this estimate to update the extracted source counts. We iterate this procedure 
until we have an estimate of the source crowding contamination for all sources in the field. 
This analysis shows that the contamination is negligible (< 1% of extracted counts) for 
all MSPs except O and R, which each have ~ 13 additional background counts due to 



1 http : / / cxc . harvard . edu /ciao/ 
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nearby sources. We include these in our estimate of their backgrounds in the subsequent 
analysis (Table [3]). Assuming these background estimates, we detect sources at each of the 
17 independent pulsar positions with > 5-cr significance. 



Throughout the analysis we assume a distance of 4.85 kpc to 47 Tuc (IGratton et al. 



20031 ). We apply an energy correction factor of 5.044 xlO -12 erg cm -2 s _1 or 1.42 xlO 34 erg s _1 
(0.5-2.0keV) per 1 HRC-S count/sec, which is the unabsorbed X-ray flux assuming A r n= 1.3x1 20 cm" 
and blackbody emissio n with a temperature of 0.178keV determined from WebPIMMS^ 
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jBogdanov et alihood ). In addition to the extraction region size, we list the total counts 



extracted, expected background counts and the time averaged luminosity in Table [3j 



3. Variability Analysis 

Prior to the timing analysis, we use the CIAO tool axbary to convert the event times 
to the solar system barycenter using the JPL DE200 solar system ephemeris, the Chandra 
orbital ephemeris and the radio/optical position of each MSP (see § |2]). 



3.1. Rotational Variability 

We calculated the rotational phase of each arri ving photon for ea ch MSP using the latest 



radio ephemerides corrected to X-ray frequencies (IFreire et al.ll2003l ; Freire et al., in prepa- 



ratio n) . The resulting phases were searched for variability with the Z 2 -test (IBuccheri et al. 



1983) . where n is the optimal number of harmonics as determined from the H-test (jde Jager et al 



19891 ). The variable Z 2 has a probability density function distributed as \ 2 with 2n degrees 



of freedom. We list the value of this variable, the detection significance in equivalent Gaus- 
sian a and n in Table [3J Only 47 Tuc D, O and R show variability above the ~ 4-cr level. 
Their folded pulse profiles are shown in Figured! All three pulsars' profiles are characterized 
by large duty cycles. 47 Tuc O shows evidence for two peaks centered at phases, ~ 0.0 and 
ps 0.4 with widths of 5(f) ~ 0.2 and S(p ~ 0.4, respectively. Only a single peak is evident 
in the profiles of 47 Tuc D and R centered at ~ 0.45 with 8(j) ~ 0.25 and <fi ~ 0.5 with 
5<p p» 0.3, respectively. 

As seen in Table [3J 5 of the 19 MSPs are detected with marginal significance, 2.8-3.3- 
o . Given the relatively large size of our MSP sample, we can quantify the significance of 
these marginal detections. We are free to choose the significance level with which we call an 



2 ht tp : / /heasarc . gsfc . nasa.gov/Tools / w3pimms . html 
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MSP 'variable'. Once we choose this level, the problem becomes one of binomial statistics 
where each MSP represents an independent measurement for variability. Choosing the 99% 
percent confidence level for Z\ (corresponding to 2.58-cr) allows us to identify 8 'variable' 
sources when only « 0.2 are expected if the MSPs were drawn from a random distribution. 
The probability of 8 of 19 trials being labeled 'variable' is 6.7xl0~ 12 . The probability of 
one or more of these being false detections is 17.4%. Since the confidence level at which we 
label an MSP as 'variable' is arbitrary, we list a range of confidence levels and the number 
of corresponding detections with the binomial and false detection probabilities in Table HI 
For the remainder of our analysis we will adopt the 99% confidence level. This results in 
47 Tuc D, E, F, H, O, Q, R and S being labeled as 'variable.' Their profiles are shown in 
Figure [H 

The pulsed fraction, f p , for each MSP can be determined using two steps. First, 
we estimate the DC ( unpulsed) level with the non-parametric bootstrapping method of 
Swanepoel et al.l (119961 ) . The advantage of this technique is that it works on the raw phases 
without the need to construct a phase histogram or know the pulse shape a priori. This level 
is shown by the solid line in Figure [1] with ±l-cr errors. Next, we correct for the fact that 
the DC component includes both unpulsed photons from the MSP and background photons. 
Note that for MSPs F and S we consider 50% of the background subtracted source counts to 
be unpulsed background photons. The expected background contribution to the pulse profile 
is denoted by the dashed line in Figured] (see § Ej). The pulsed fraction determined by the 
bootstrapping method can be related to the true pulsed fraction by f p = f P) haotN t J (N t — N b ), 
where / Pi boot is the pulsed fraction determined by the bootstrapping method, N t is the total 
number of extracted counts and iV b is the total number of background counts contributing 
to N t . Both N t and Nj, can be found in Table [31 The pulsed fractions derived in this manner 
are also listed in Table [31 

We calculate upper limits on the pulsed fractions of the remaining 11 MSPs with tim- 
ing solutions (note that W has only a preliminary solution) with Monte Carlo simulations 
assuming two scenarios. First, we make the conservative assumption that the underlying 
pulse shape is sinusoidal. For each MSP and pulsed fraction we simulate 500 light curves 
with a total number of counts N t = Nj, + N s , where N s is the number of source counts. If 
the pulsed fraction is f P)S i n e, then N s consists of (1 — f PtS i ne )N s unpulsed counts and / P]Silie Af s 
pulsed counts. For the two unresolved MSPs that do not show pulsations (G and I) we 
assume that 50% of the source photons constitute an unpulsed background. We determine 
the 90% confidence upper limit on the pulsed fraction as the value of / P]Sijie at which 450 
synthetic light curves have values of Z\ ^ 27.4 (corresponding to 5-er). The result is that 
only J, L, W and Y have sufficient counts to be constrained in this way. 
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Motivated by the close correspondence of the radio and X-ray pulse profiles of magne- 
tospherically dominated MSPs (see § HJ, we determine a second set of upper limits based 
on the assumption that the underlying X-ray pulse has the same morphology as the radio 
pulse. We denote this pulsed fraction as / Pira dio- In the case of the weak radio pulsars, N, R, 
T, W and Y it was necessary to model the radio pulse(s) with Gaussians and use these as 
the assumed X-ray pulse shape. The radio pulse profiles for remaining MSPs have sufficient 
signal-to-noise to be used directly. We then perform the Monte Carlo simulations with the 
technique described above to determine the 90% confidence upper limits. The results can be 
found in Table [31 It is possible to constrain the X-ray emission with the same morphology 
as the radio pulse from 7 pulsars in this way. 



3.2. Orbital Variability 

The orbital periods of the 12 binary 47 Tuc MSPs span the range 0.07-2.4 days. In 
order to search for variability during these orbits we calculated the orbital phase of each 
arriving photon for each MSP and constructed histograms with 5, 10 and 20 bins. We then 
corrected each bin for the variation of exposure time during that particular orbital phase 
so that we have a histogram of counts per second per bin and subtracted the expected 
background contribution. We search for variability by computing \ 2 between the histogram 
and a constant count rate. We found that none of the 47 Tuc binary MSPs showed significant 
orbital variability. 

Substantial X-ray eclipses characterized by a complete disappearance of hard (> 2keV) 
photons and a decline in soft (< 2 keV ) photons for « 30% of the orbit have been reported 



for 47 Tuc W (IBogdanov et al.l 120051 ) . Thus, the lack orbital variability in 47 Tuc W is 
surprising. We reduced archival ACIS-S data of 47 Tuc W from ObsIDs 2735, 2736, 2737, 
and 2739 in order to quantify the properties of the eclipse in the HRC-S bandpass. We 
reprocessed the level 1 event files to make use of the latest calibration and filtered periods of 
high background flaring (i.e. periods with count rates > 3-cr above the mean full- frame rate). 
This resulted in 251.3 ksec of exposure time. We computed the orbital phase of each photon 
within a 1" radius of the position of 47 Tuc W in the energy range 0.3-2.0 keV with the 
preliminary timing solution (see § |2J). The X-ray eclipse is evident in the histogram shown 
in Figure [2h at phase, ~ 0.2, where <fi = at the time of the ascendin g node. This X-ray 



eclipse timing agrees well with observed radio eclipse from = 0.1-0.4 (jCamilo et al.l 12000 



Edmonds et al.l 120021 ) . The variability is significant at the 4.3-cx level from the measured 
value of Zf = 30.1 and we determine an orbital modulation of f p = 36 ± 9% using the 
non-parametric bootstrap method (see § 13.11) . This is consistent with the 90% confidence 
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limit of / P)S ine < 48% on a sinusoidal signal in the HRC-S time series of 47 Tuc W (see § 13. ip . 



3.3. Aperiodic Variability 



To search for aperiodi c varia bility in the 47 Tuc MSPs, we have applied the Bayesian 
blocks algorit hm of IScarrie! Jl998h as im plemented by the Interactive Spectral Interpretation 
System (ISIS; Houck fc Denicola] 2000 ). The algorithm determines the optimal decomposi- 
tion of the light curve into constant count rate segments based on a parametric maximum 
likelihood model of a Poisson process. The raw (unbinned) events are divided into 'blocks' 
and the odds ratio that the count rate has varied is computed. If variability is found, each 
'block' is further subdivided to characterize the structure of the light curve (e.g. step-function 
variation, flaring, etc.). We could not identify any intra-observation variability from 17 MSPs 
using an odds ratio corresponding to 68% chance that any variability is real. In addition, the 
inter-observation count rates for each MSP derived from this process do not show significant 
variability over the ~ 20 day span of our observations. Thus, the X-ray emission from the 
47 Tuc MSPs is stable on timescales ranging from minutes to days. 



3.4. ACIS-S vs. HRC-S Comparison 



With such a large sample of constant luminosity X-ray sources, we can compare the count 
rates between the ACIS-S and HRC-S for soft thermal sources. In Figure [3] we compare the 
count rate for the two detectors using each of the 17 independent MSP detections. For the 
ACIS-S count rate, we summed the values in the 0.3-0.8 keV and 0.8-2.0 keV bands listed 
by Heinke et al. (2005). The relation between count rates in the two different detectors was 
/(HRC - S) = (0. 43 ± 0.024) x /(ACTS - S). This is consistent with the conversion from 
HRC-S to ACIS-I feutledge et aPboO^ . 



4. HRC-S Timing Accuracy 



Th e accuracy of the HRC -S time tags was demonstrated to be ±12 /xs in an observation 
of M28 ORutledge et al.ll2004j ). We have investigated several issues to ensure that the HRC-S 
has sufficient accuracy to detect MSPs. In addition to accounting for a leap second that 
occurred during the middle of the observing span, we examined the effect of telemetry 
saturation, different solar system ephemerides and analyzed a recent calibration observation 
of the globular cluster M28. 
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4.1. Telemetry Saturation 

The maximum telemetered full-field, unfiltered count rate for the HRC-S is 184 counts 
.-'(Chandra Proposers' Observatory Guideft. At rates above this, a decreasing fraction of 
all counts will be telemetered. The data will be affected if this count rate is exceed during 
the 2.05 s full-frame readout time. The effect of telemetry saturation on timing certainty 
is that, when telemetr y is saturated, not a ll events are telemetered back to Earth. Due to 



the HRC wiring error (ITennant et al.ll200ll ). the Nth event detected by the HRC-S has its 
time assigned to the iV+lst event 8 . The true time series can be reconstructed if both are 
telemetered; if either is not, then the wrong time will be assigned to one event (either the 
iVth event, or the iV — 1st event). 

To test the extent to which this saturation impacts our observations we constructed a 
light curve of our entire dataset binned in 2.05 s intervals. We found that only « 0.1% of the 
bins exceeded the maximum count rate. Thus telemetry saturation will effect only fh 0.1% 
of the counts in any given MSP, which is negligible. 



4.2. Ephemeris Comparison 

We performed a preliminary extraction of data using DE405 and compared the timing 
precision of data corrected to DE200 using one of the brightest sources in the field. We 
found the entire observational period was offset by £de405 — ^DE2oo = 1.3809 ms at the start 
of the observational period, decreasing monotonically to 1.3719 ms at the end of the obser- 
vation period. Thus, there was an average direct offset between the photon time-of-arrivals 
(TOAs) in the two ephemerides of «1. 3764 ms, and a range of variation of ~ 9 /xs. This 9 /is 
therefore amounts to the relative timing uncertainty due to the adopted ephemeris, compa- 
rable to the uncertain ties in HRC digitization (±5 (is) and Chandra clock stability (±5 /is) 



flRutledee et all 12004 ). 



4.3. s/c Clock Stability 

A calibration observation of PSR B1821— 24A was performed on May 27, 2006 starting 
at 12:30 UT for 40887 sec in order to evaluate the stability of the HRC-S clock. A complete 
analysis is beyond the scope of this work. However, a search of the data using the known radio 



3 http : / / cxc . harvard . edu /proposer/POG /index, html 
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timing ephemeris (IRutledge et al.1 120041 ) clearly detects the 3.05 ms pulsar with Z\ = 330 
corresponding to a detection significance of ~ 18-c. Thus, we conclude the HRC-S clock has 
remained sufficiently stable to detect MSPs. 



Discussion and Conclusions 



Chandra HRC-S observations of 47 Tuc have allowed us to study a relatively large 
sample of 19 MSPs on millisecond to week timescales. We find that the MSPs in 47 Tuc 
uniformly show very low levels of variability on all scales probed. We have sufficient statistics 
to meaningfully constrain (under the assumption that their X-ray pulse profiles match the 
radio) or measure the rotational modulation of 15 MSPs. Eight of these objects have low 
pulsed fractions, f p ^ 50%. MSPs 47 Tuc D, O, and R each have pulsations detected at ^ 
4-er significance with relatively large pulsed fractions, f p ^ 60%, which are similar to the 
levels seen from luminous MSPs dominated by non-thermal emission (e.g. PSRs B1821— 24A 
and B1937+21). However, the pulse profiles of these objects (see Figure [T]) are characterized 
by large duty-cycle features tha t do not resemble the sharp , low duty-cycle profiles seen in 



the non-thermal MSPs (e.g. see lBecker fc Aschenbachll2002l ). 



The existing ACIS data show that the 47 Tuc MSPs ha ve fairly homogeneous spectro- 
scopic properties (IGrindlay et al.ll2002l ; iBogdanov et al.ll2006l ). All but 3 of the 47 Tuc MSPs 
are characterized by l-3xl0 6 K thermal spectra with low luminosities in a narrow range, 
Lx ~ 10 30-31 erg s -1 , and have small emission radii, i? e ff ~ 0.1-3.0 km. The low level of mea- 
sured variability presented here indicates that rotational averaging do es little to affect these 
value s, which agree with the predictions of polar cap heating scenarios (IHarding fc Muslimov 
20021 ). Thus we conclude that, unlike the luminous non-thermal MSPs, the vast majority 
of the X-ray emission from the 47 Tuc MSPs is created by the heating of the neutr on star 
polar cap by a return current of r e lativistic particles produced in the magnetosphere (lArons 



198ll : IHarding fc Muslimov! l200ll . 12002 ; iGrindlav et al.l 12002 ; IBogdanov et al. 



20061 ). For 



older MSPs with very short spin periods and low magnetic fields, like those in 47 Tuc, the 
main source of the e ± -pair production is thought to be through inverse Compton scattering 
of thermal X-rays from the neutron star surface off of electrons in the pulsar magnetosphere 
(IHarding k Muslimovi [20021 ) . 



Only the radio-eclipsing binaries 47 Tuc J, O and W show power-law spectral com- 
ponents that contribut e significantly (70%, 50% and 75%, respectively) to their total flux 
(IBogdanov et al.l 120061 ) . The lack of strongly pulsed emission in 47 Tuc J and W suggests 
that the X-ray emission does not arise in the neutron star magnetosphere, but instead is 
likely the consequence of an intra-binary shock. This is in agreement with conclusions based 
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on orbital phase resolved spectroscopy of 47 Tuc W by iBogdanov et al.l (120051 ). Conversely, 
the current data do not conclusively identify the origin of X-rays from 47 Tuc O, which shows 
significant pulsations. The X-rays from an intra-binary shock are not expected to be modu- 
lated at the rotational period of the MSP, so the measured pulsed fraction, f p = 83 ±21%, is 
only marginally consistent with the 50% spectroscopic allocation of X-rays due to a shock. 
In addition, its large duty-cycle does not immediately imply that polar cap heating is the 
source of the pulsed X-rays, since broadl y beamed magnetospheri c emission viewed off-axis 



would appear to have a large duty-cycle (IBecker fc Triimper 



1999|). 



The apparent non-detection of low duty-cycle pulsars is significant in comparison with 
the pulse profile of B 182 1—24 A. If all 19 MSPs had X-ray pulse profiles identical to that of 
PSR B1821— 24A, all would have been detected with ^ 7-a significance (which we find for the 
lowest SNR MSP, 47 Tuc T). Those MSPs with higher count rates would have been detected 
with even greater significance. The implication is that PSR B1821— 24A has an unusually 
low duty cycle for a MSP. If we assume that low duty-cycle MSPs make up a fraction / 
of the globular cluster population, then the non-detection of even 1 such pulsar in 47 Tuc 
implies that B1821-24A-like pulsars comprise / <20% of the MSP population in GCs (90% 
confidence limit). This limit could be much lower, if the intrinsic distribution of duty cycles 
in magnetospheric pulsars is lower than that of B1821— 24A, for example; however, there 
seems to be little observational work quantifying the distribution of duty cycles of observed 
pulsars. 
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Table 1. X-ray Properties of Millisecond Pulsars Outside of 47 Tuc. 



PSR 


P 


d 


T 


logE 


log L x 


f P 


Refs. 




(ms) 


(kpc) 


(Gyr) 


(ergs -1 ) 


(ergs -1 ) 


(%) 






B1937+21 


1.56 


3.57 


0.23 


36.04 


33.15 


54 ± 7 


1,2 


B1957+20 


1.61 


2.49 


2.22 


35.04 


31.81 


< 60 


1,3 


J0218+4232 


2.32 


2.67 


0.48 


35.38 


32.54 


59 ± 7 


4,5 


B1821-24A (M28) 


3.05 


5.5 a 


0.03 


36.35 


33.22 


85 ± 3 


6,7 


J0751+1807 


3.48 


1.15 


7.08 


33.86 


30.84 


52 ± 8 b 


8,9 


J0030+0451 


4.87 


0.32 


7.71 


33.53 


30.40 


69 ± 18 


10,11 


J2124-3358 


4.93 


0.25 


6.01 


33.83 


30.23 


56 ± 14 


1,12 


J1012+5307 


5.26 


0.84 


4.86 


33.68 


30.38 


77 ± 13 b 


13,9 


J0437-4715 


5.76 


0.14 a 


4.91 


33.58 


30.47 


40 ± 2 


14,12 


J1024-0719 


5.16 


0.39 


27 


32.93 


29.30 


52 ± 22 


1,12 


J1744-1134 


4.07 


0.36 a 


9.1 


33.62 


29.49 c 




1,3 


J0034-0534 


1.88 


0.54 


6.03 


34.48 


29.60 c 




12 


No High Time Resolution Imaging 


B1620-26 (M4) 


11.08 


1.73 a 


0.26 


34.28 


30.08° 




15,16 


J1740-5340 (NGC 6397) 


3.65 


2.55 a 


0.34 


35.15 


30.9 C 




17,18 


J1911-6000C (NGC 6752) 


5.28 


4.1 a 


38.1 


32.77 


30.34 c 




19 


J2140-2310A (M30) 


11.02 


9.0 a 


>0.08 


<34.79 


30.64 c 




20 



a Accurate distance measurement. 

b Detection significance is low. 

c X-ray luminosity in the 0.5—2.5 keV band. 



Note. — All distances are estimated from th e pulsar dispersion measures and the model of Galactic dis- 
tribution of free electrons llCordes fc Lazicll2002h, e xcept where noted. X-ray luminosities are quoted in the 
0.2-10 keV band as adopted from Table 1 of IZavlinl l|2006h and references therein, except where noted. Pulsed 
fractions are quoted in roughly in the HRC-S band (0.1-2.0 keV), but see references for the specific bandpass. 
The spectra of MSPs above the first horizontal line are dominated by non-thermal X-ray emission. Those 
below the line have significant thermal components or are indeterminate (and thu s presumed to be t h ermal ) 



in nature. References: 1 — iToscano et al 
4 - iNavarro et all lll995l) , 5 



Webb et al 



_[RjrUedge^^alT 

Nicci'l id.l |2005h , 9 - Webb et al.l j2004b f)7 10 - Lommen et all (120001 ) . 11 -iBecker fc Aschenbacnl 



|NicastroetalJ 
^BeckeT^TldT 



- Becker 



Triim 
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Zavlin 
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20041) , 17 - iD'Amico et al.1 ibOQlO s - iGrindlav et al.l d2002h . 
20041) . 
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Table 2. HRC-S Observations. 



ObsID Obs. Start Obs. Start Exp. Time 
(UT) (MJD) (sec) 



5542 


2005 


Dec 


19.29 


53723. 


.79 


51918.2 


5543 


2005 


Dec 


20.62 


53725. 


12 


53962.6 


5544 


2005 


Dec 


21.98 


53726. 


.48 


52036.6 


5545 


2005 


Dec 


23.21 


53727. 


71 


54203.5 


6237 


2005 


Dec 


24.59 


53729. 


09 


51920.8 


6238 


2005 


Dec 


25.88 


53730. 


.38 


50887.6 


5546 


2005 


Dec 


27.23 


53731. 


73 


51939.3 


6230 


2005 


Dec 


28.57 


53733. 


.07 


52401.0 


6231 


2005 


Dec 


29.91 


53734. 


.41 


48963.7 


6232 


2005 


Dec 


31.22 


53735. 


.72 


49139.0 


6233 


2005 


Jan 


2.24 


53737. 


.74 


103433.2 


6235 


2005 


Jan 


4.17 


53739. 


.67 


51932.1 


6236 


2005 


Jan 


5.48 


53740. 


.98 


54729.8 


6239 


2005 


Jan 


6.92 


53742. 


.42 


52241.6 


6240 


2005 


Jan 


8.10 


53743. 


60 


54178.4 



Table 3. HRC-S Derived X-ray Properties of the 47 Tuc Millisecond Pulsars. 



Name 


Total a 


Background 


Extraction 


log(L x ) b 


Harmonics c 




Significance 


fp, radio 


' fp, sine 






Counts 


Counts 


Radius 


(0.5-2.0 keV) 


(n) 












AT Tii,- P 
4 ( 1UC \_j 


"1 7Q 


y4.o 


1. o 


on in 
OU. 1Z 


1 

1 


9 K 
Z.O 


l.o 




, 1 nn 




AT Tii/> T\ 
4 ( 1UC D 


9 "1 


y4.o 


1. o 




n 
Z 


91 7 
Zl. f 


Q ("1 






ID ± Zl 


AT Tur T?S 
4 ( 1UC .Cj 


ZD4 


y4.o 


1. o 


on a o 
OU. 4o 


1 

1 


1Z. 


o . o 






OU it ly 


/17 Tii,-. T e 
4 ( 1UC r 


-1 1 Q 
41o 


on 7 
oU. i 


1"9 
1 . z 


^-"3n 7c; 


1 

1 


in 1 


9 ("1 






Z It IZ 


47 Tnr> pf 
4/ _LUC \j 


'^99 


sn 7 

OU. ( 


1"9 
1 . z 


xOU.Ol 


1 
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4/ 1UC ri & 
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y4.o 


1. o 
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OU. 14 


n 
Z 


1 O 7 
1Z. / 


Z.o 






oft _l_ on 
ZD it ZU 


/IV Tn/- Ti\<7 
4/ 1UC 1 ,a 


Q99 


on 7 
oU. i 


i//n 
1 . Z 


<loU.Dl 


1 


o.u 


9 Q 
Z.O 


<ol 


^- -i nn 




47 Tuc js 


266 


94.8 


1"3 


30.46 


1 


6.5 


2.3 


<38 


<77 




47 Tuc L 


342 


35.9 


c/s 


30.71 


2 


10.0 


2.4 


<50 


< 49 




47 Tuc M 


151 


94.8 


1'.'3 


29.98 


1 


2.2 


1.4 


< 100 


< 100 




47 Tuc N 


186 


94.8 


1'.'3 


30.19 


1 


4.8 


2.0 


<73 


< 100 




47 Tuc O g 


431 


94. l h 


1"2 


30.77 


3 


39.1 


5.1 






81 ±21 


47 Tuc Qs 


186 


94.8 


1'.'3 


30.19 


6 


28.3 


3.1 






83 ±42 


47 Tuc RS 


288 


80.4 h 


1"1 


30.57 


2 


24.1 


4.1 






63 ±29 


47 Tuc S C 'S 


413 


80.7 


1"2 


<30.75 


8 


33.4 


3.1 






20 ±15 


47 Tuc TS 


133 


80.7 


1"2 


29.94 


1 


1.2 


1.1 


< 100 


< 100 




47 Tuc U g 


193 


94.8 


1"3 


30.22 


1 


0.1 


0.7 


< 100 


< 100 




47 Tuc Ws 


433 


80.7 


1"2 


30.77 


1 


0.8 


1.0 


<48 


<48 




47 Tuc YS 


218 


94.8 


1"3 


30.32 


1 


(5.5 


2.3 


<40 


<96 





-J 



a This number includes the expect background counts listed in the subsequent column. 
b X-ray luminosity (logarithm of erg s _1 ) derived in the band 0.5— 2.0keV. 
c See the text (§ I3,lj l for a description of these parameters. 

d <100% means the undetected pulsation is consistent with 100% pulsed signal, and therefore is unconstrained by simulations. 

c 47 Tuc F and S have overlapping positions. The total counts represent all photons extracted from the l'/2 extraction radius 
and the background counts are only those expected from a uniform background. 

f 47 Tuc G and I have overlapping positions. The total counts represent all photons extracted from the l'/2 extraction radius 
and the background counts are only those expected from a uniform background. 

^Binary MSP. 

h Includes an estimate of the contamination due to source crowding (see § 13,11 . 
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Table 4. Summary of 47 Tuc MSP Detection Significance. 



Pulsation Detection 


Number of 


Binomial 


False Detection 


Threshold 


Detections 


Probability 


Probability 


(1) 


(2) 


(3) 


(4) 


99% 


8 


6.9xlfr 12 


17.5% 


99.73% 


6 


f.Oxfcr 11 


5% 


99.947% 


3 


i.4xicr 7 


1% 



Note. - (1) - Confidence level of Z\ above which a 47 Tuc 
MSP is labeled 'variable.' (2) - The number of 'variable' MSPs for 
the given confidence level. (3) - The binomial probability that the 
number of 'variable' MSPs is due to chance. (4) - The probability 
that one of the 'variable' MSPs is a statistical fluctuation. 
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Fig. 1. — Pulse profiles of the variable 47 Tuc MSPs with a typical error bar (two periods 
are shown for clarity). The histograms were constructed with 20 bins per period. The solid 
horizontal line denotes the DC (unpulsed) contribution to the pulse profile as determined by 
the nonparametric bootstrapping algorithm with the associated 1-a uncertainties denoted 
with dashed-dotted lines. The dashed line denotes the estimated contribution to this level 
due to the uniform background and source crowding (50% of the background subtracted 
counts are attributed to this total for MSPs F and S; see § 13. II for details). 
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Fig. 2.— Orbital Profiles of 47 Tuc W from ACIS-S (upper) and HRC-S (lower). The error 
bars are 1-a and histograms contain 20 bins per period. 
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Fig. 3. — HRC-S/ACIS-S count rate comparison using the 17 independent 47 Tuc MSP 
detections. The line represents the best linear fit to the data. 



